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要 : 干旱 半 干 旱地 区 急需 高 分 状 率 的 土壤 盐 度 图 用 于 显示 盐 度 空间 分 布 的 细微 变化 ,指导 盐 渍 


化 区 域 和 潜在 盐 涡 化 区 域 制 定 土地 资源 管理 政策 和 水 资源 管理 政策 ,防止 土壤 进一步 退化 ,保障 
农业 经 济 可 持续 发 展 和 粮食 安全 生产 。 基 于 PlanetScope 影像 ,提取 植被 光谱 指数 和 土壤 盐 度 指 


数 ， 


共计 21 个 变量 ,将 其 输入 装 袋 回归 (Bootstrap aggregating,Bagging) 算 法 中 ,构建 了 土壤 盐 度 预测 


模型 Model-I ;使 用 最 相关 最 小 元 余 (Max-relevance and min-redundancy,mRMR ) 方 法 第 选 特征 变 


量 
E, 


将 其 输入 Bagging 中 ,构建 了 土壤 盐 度 预测 模型 Model-I ,使 用 野外 采样 数据 来 辅助 建 模 并 进行 


验证 。 通 过 模型 评价 指标 对 Model-I 和 Model-I 进行 评 估 。 结 果 表 明 :Model-IT 的 预测 性 能 优 
Model-I (验证 集 决 定 系 数 为 0.66, 均 方 根 误差 为 18.00 dS «m! ,四 分 位 数 的 相对 预测 误差 为 3.21 ) 
mRMR 有 效 降低 了 多 维特 征 宛 余 问 题 。PlanetScope 影像 结合 mRMR 方法 成 功 绘制 了 高 分 辩 率 了 


> 


p Eo 


度 图 ,提供 了 更 详细 的 土壤 盐 度 空间 分 布 信 息 , 研 究 结果 对 利用 PlanetScope 数 据 监测 土壤 盐 


渍 化 信息 起 推动 作用 。 
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土壤 盐 省 化 阻碍 了 植被 对 水 分 的 吸收 ,加 速 了 
土壤 表面 有 机 质 分 解 ,降低 了 土壤 肥力 和 土壤 团 只 


化 ; 数字 制图 技术 


的 有 效 指标 ,但 在 水 平和 垂直 方向 呈现 出 快速 的 时 
空 动态 变化 ”。 因 此 ,通过 制作 十 壤 盐 度 分 布 图 能 


体 稳定 性 ,是 土壤 退化 的 主要 因素 之 一 小。 全 球 范 
围 内 有 7.7%~11.5% 的 灌溉 土地 受到 盐 渍 化 威胁 ,这 
意味 着 过 量 的 可 溶性 盐分 累积 ,对 作物 生产 和 农业 
可 持续 发 展 带 来 严重 影响 小 。 我 国 发 生 盐 渍 化 的 地 
区 分 布 广泛 ,其 中 新 疆 地 区 的 盐 渍 十 类 型 最 多 , 盐 
渍 土 总 面积 达 847.6x10: hm , 占 新 疆 地 区 现 有 耕地 
面积 的 31.1%”。 充 足 的 耕地 资源 是 保障 粮食 生产 
和 安全 的 基石 ,但 不 断 加 剧 的 土壤 盐 渍 化 每 年 将 
1%~2% 的 耕地 转化 为 盐 渍 地 ,给 粮食 生产 带 来 了 严 
重 限制 和 挑战 所 。 土 壤 盐 度 是 评估 盐 涡 化 严重 程度 
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帮助 决策 者 了 解 土壤 盐 度 存在 的 巨大 时 空 差 异 ,为 
潜在 盐 渍 化 区 域 制定 优化 的 土地 资源 管理 和 水 资 
源 管理 政策 ,为 盐 涡 化 区 域 制 定 优化 的 土壤 改良 政 
策 , 从 而 应 对 新 疆 地 区 土地 持续 性 退化 难题 。 

20 世 纪 70 年 代 , 新 疆 地 区 通过 第 二 次 全 国土 地 
调查 获取 了 土壤 盐 度 数据 ,使 用 矢量 空间 数据 结构 
记录 其 空间 信息 ,这 种 数据 记录 方式 难以 准确 地 反 
映 土壤 属性 的 空间 异 质 性 ,也 难以 捕获 盐 度 信息 的 
动态 变化 ,无 法 满足 新 疆 地 区 对 土壤 盐 涡 化 管理 和 
预防 的 需求 "。 研 究 证 明 ,光谱 学 知识 结合 遥感 技 
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术 可 以 快速 昌 定 量 地 获取 土壤 盐 度 信息 ,并 成 功用 
于 绘制 高 分 辨 率 的 土壤 盐 度 图 ,提供 了 更 全 面 的 土 
坏 盐 渍 化 空间 信息 ”。 盐 渍 化 土壤 相 较 于 非 盐 渍 化 
土壤 ,在 可 见 光 到 短波 红外 波长 范围 内 具有 更 高 的 
光谱 反射 率 , 随 着 土壤 盐 度 增加 ,光谱 反射 率 进 一 
步 提高 ,在 近 红 外 波长 附近 出 现 具有 诊断 性 的 吸收 
M" 此 , 归 一 化 盐 度 指数 (Normalized differ- 
ence salinity index, NDSI) 和 多 种 盐 度 指 数 (Salinity 
index, ST) 被 广泛 用 于 评估 土壤 盐 度 。 植 物 对 土壤 
盐 碱 变化 的 耐 受 ,通过 生长 状况 进行 响应 ,长 势 不 
好 的 植被 光合 作用 减弱 ,叶绿素 含量 变 少 , 可见光 
范围 内 反射 率 增强 , 近 红 外 范围 内 反射 率 降 低 ”。 
此 ,上 归 一 化 植被 指数 (Normalized differential vege- 
tation index,NDVI) 和 土壤 调节 植被 指数 (Soil adjust- 
ed vegetation index,SAVI) 等 多 种 植被 指数 (Vegeta- 
tion indices, VIs) 被 广泛 用 于 评估 土壤 盐 度 的 间接 
指标 。Gorji 等 ' 引 基于 Landsat8 OLI 和 Sentinel-2A 
MSI 影 像 ,成 功 提取 了 与 土壤 电导 率 相 关 的 SL, 并 结 
合 多 元 线性 回归 技术 成 功 绘制 了 乌 尔 米 耶 湖 流 域 
土壤 盐 度 图 。Fernandez-Buces 等 "基于 Lansat 
ETM 机 载 近 红外 和 实测 光谱 数据 获取 了 不 同 植被 
和 裸 地 的 反射 率 , 提 取 并 调整 了 NDVL, 并 计算 了 纤 
合 光 谱 指 数 (Combined spectral response index , COS- 
RI) ,以 制作 出 高 精度 的 原 墨 西 哥 特 斯 科 科 湖区 域 
HRERL Ding 等 ”使 用 近 端 遥 感 装置 EM38 测 量 表 
观 电导 率 数 据 ,基于 Landsat TM 影像 提取 了 SI, 
NDSI 亮度 指数 (Britness index, BI) .NDVI 和 差 值 植 
被 指数 (Difference vegetation index , DVI) ,将 表 观 电 
导 率 数据 和 遥感 指数 信息 应 用 于 回归 克 里 金 插值 
法 ,成 功 评估 了 土壤 盐 度 并 获得 了 满意 的 结 

20 RR ,机 器 学 习 算 法 被 广泛 使 用 ,在 土壤 


度 预 测 模型 时 ,特征 变量 的 选择 至 关 重 要 。Allbed 
等 ”2” 的 研究 结果 表明 植被 指数 和 盐分 指数 等 特 
征 变量 可 以 有 效用 于 土壤 盐 渍 化 监测 研究 ,但 是 这 
些 特征 变量 均 是 通过 波段 间 的 运算 得 到 ,造成 了 特 
征 变 量 间 不 同 程度 的 信息 匈 余 。 最 相关 最 小 宛 余 
(Max-relevance and min-redundancy , mRMR ) 特 征 选 
择 算法 具有 较 低 的 复杂 度 , 与 当前 常用 的 特征 选择 
算法 相 比 具有 原理 简单 和 对 数据 要 求 低 的 优点 ， 
mRMR 具有 的 优势 提升 了 土壤 盐 度 预测 模型 的 预测 
效率 2 。Peng 等 2 使 用 Pearson 相关 分 析 筛 选 相 关 
性 高 的 特征 变量 进行 土壤 盐 渍 化 监测 ,并 没有 考虑 
选择 特征 间 存 在 的 信息 元 余 。 马 慧 琴 等 ”在 筛选 
特征 用 于 构建 遥感 监测 模型 时 ,对 比 了 mRMR 与 常 
用 相关 分 析 方 法 ,人 研究 结果 表明 mRMR 方法 比 常 用 
相关 分 析 方 法 更 具 优 势 。 

目前 ,遥感 影像 的 中 低空 间 分 辩 率 仍然 是 影响 
土壤 盐 度 图 精度 的 主要 因素 ”。 当 受 盐 影响 区 域 
小 于 像素 大 小 时 ,会 限制 土壤 盐 度 图 对 土壤 盐 度 进 
行 准确 评估 的 能 力 ,空间 分 辨 率 低 于 5 m 的 多 光谱 
传感器 正成 为 绘制 土壤 盐 度 图 不 可 或 缺 的 装置 ” 。 
Bouaziz 等 ”证实 了 基于 MODIS terra 数据 提取 的 光 
谱 指 数 与 实测 数据 能 有 效 绘制 土壤 盐 度 图 ,但 较 低 
的 空间 分 辩 率 是 造成 二 者 相关 性 较 差 的 主要 原 
o Eldeiry 45 ^" FE $$ T IKONOS 与 Landsat 影像 在 
农业 地 区 评估 土壤 盐 度 的 性 能 ,结果 表明 高 空间 分 
辨 率 能 获取 更 准确 的 评估 结果 。 自 2016 年 以 来 ， 
Planet 公司 运营 的 地 球 成 像 小 型 卫星 星座 Planet- 
Scope ,能 够 以 较 高 的 空间 分 辨 率 (<5 m) 和 时 间 分 辩 
率 (1~2 d) 获 取 影 像 数据 。 在 2019 年 4 月 Planer- 
Scope 携 带 了 第 三 代 传 感 器 PSB.SD 能 获取 可 见 光 近 
红外 范围 内 的 8 个 波段 YW。PlanetScope 获 取 的 影像 


科学 领域 ,基于 机 器 学 习 算 法 绘制 了 多 种 数字 土壤 
地 图 ,如 Farifteh 等 "基于 随机 森林 方法 (Random 
forest, REF) 开发 了 绘制 土壤 有 机 碳 浓度 和 储量 数字 
图 的 框架 。Ma 等 基于 RF 和 极端 梯度 提升 等 机 
带 学 习 方 法 绘制 了 渭 库 绿洲 土壤 盐 度 图 ,有 效 解 决 


与 Landsat8 OLI ^" 和 Sentinel-2A MSI ^" fH EG , HA 
更 高 的 空间 分 辨 紊 和 时 间 分 辨 紊 ,而 且 国 内 利用 
PlanetScope 小 型 星座 卫星 影像 进行 土壤 盐 度 制图 
的 研究 并 未 广泛 展开 。 促 进 高 空间 和 高 时 间 分 辩 
率 的 遥感 影像 用 于 土壤 盐 度 监测 研究 ,能 获取 更 翔 


了 土壤 和 环境 因子 间 的 非 线性 问题 。Wang 等 ”在 
绘制 不 同 区 域 的 土壤 盐 度 图 时 ,比较 了 5 种 机 器 学 
习 算法 ,发 现 基于 机 天 学 习 算 法 进行 土壤 盐 度 预测 
研究 时 需要 用 户 定义 的 参数 更 少 ,拥有 更 高 的 计算 
效率 ,能 够 处 理 数值 序数 和 离散 预测 变量 等 优点 。 
然而 ,在 构建 高 效 、 准 确 和 泛 化 程度 高 的 土壤 盐 


实 的 空间 信息 动态 变化 的 盐 度 信息 和 更 精确 的 评 
佑 结 

本 文 以 渭 干 河 - 库 车 河流 域 三 角 洲 绿洲 (简称 
渭 库 绿洲 ) 为 研究 区 ,基于 PlanetScope 卫星 的 原始 
波段 数据 和 波段 运算 所 得 的 光谱 指数 , 共 21 个 特征 
变量 。 采 用 mRMR 方 法 筛选 重要 特征 作为 自 变量 
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输入 Bagging 回归 算法 ,建立 了 土壤 盐 度 预测 模型 ， 
成 功 绘制 了 研究 区 高 分 辨 率 土壤 盐 度 分 布 图 。 这 
为 新 疆 地 区 制定 管理 改善 土壤 盐 碱 策 略 ,维持 绿洲 
可 持续 发 展 和 保障 粮食 安全 提供 重要 参考 。 


1 材料 与 方法 


1.1 研究 区 概况 

渭 库 绿洲 ,属于 南 疆 绿洲 区 。 地 势 北 高 南 低 、 
地 貌 复 杂 多 样 ,气候 干燥 、 节 发 强烈 降水 稀少 ,属于 
典型 的 温带 大 陆 性 干旱 气候 ”。 复 杂 的 地 形 地 貌 UR 
端的 气候 条 件 .农业 上 排水 系统 不 健全 灌溉 技术 
粗放 等 问题 ” ,导致 研究 区 存在 不 同 程度 的 盐 渍 
化 。 绿 洲 外 围 植 被 稀 玻 ,生态 系统 脆弱 ,主要 包括 
芦苇 (Phragmites australis) 、 树 柳 (Tamarix chinen- 
sis) W£ (Suaeda glauca ) MIKSERI] (Alhagi sparsifo- 
lia) 等 植物 。 实 际 人 研究 区 位 于 库 车 市 绿洲 外 围 农田 
荒漠 交错 带 , 地 理 位 置 介 于 83?23'~83?29'FE.41?23'~ 
41°27'N 之 间 ( 图 1)。 
1.2 样品 采集 与 处 理 

实验 团队 于 2021 年 6 月 23 一 24 日 在 研究 区 进 
行 了 实地 调查 ,结合 以 往 的 实地 调查 成 果 和 采样 经 


(a) 渭 库 绿洲 À 
一 研究 区 
(b) 农田 荒漠 0 10km 
采样 点 分 布 L 
À 
图 例 
o 采样 点 
T. qm 
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验 , 在 农田 荒漠 交错 带 选 择 了 84 个 具有 代表 性 的 采 
样 点 。 用 便携 式 GPS 记 录 采 样 点 的 地 理 位 置 ,以 采 
样 点 为 中 心 ,在 距离 采样 点 1 m 的 范围 内 ,以 0~10 
cm 的 深度 使 用 土壤 取样 器 均匀 的 采集 5 个 样品 , 然 
后 将 5 个 样品 混合 在 一 起 , 取 混 合 后 的 样品 约 500 
g, 装 人 有 标签 的 防水 密封 袋 。 在 实验 室 里 ,所 有 的 
土壤 样品 均 经 过 风干 .人 研磨 和 过 篇 (2.0 mm ) 等 一 系 
列 处 理 。 取 20 g 处 理 过 的 土 样 和 100 mL 纯 水 充分 
振荡 , 静 置 24 h。 在 25 %C 的 温度 下 ,用 多 参数 测量 仪 
(WTWinoLab(B) Multi3420 set B, WTW GmbH , Germa- 
ny) 测 量 静 置 后 的 浸出 液 ,得 到 电导 率 (EC) 值 。 由 
于 土壤 中 的 盐分 涂 解 在 浸出 液 中 ,浸出 液 的 含 盐 量 
与 25 SC 下 电极 的 导电 能 力 呈 正 相 关 关 系 ,通过 测量 
土壤 浸出 液 的 EC 值 , 可 以 作为 评估 土壤 盐 度 的 参 
考 指标 |。 
1.3 遥感 影像 预 处 理 与 光谱 指数 提取 

PlanetScope 是 通过 多 组 发 射 的 单独 卫星 组 成 
的 卫星 星座 , 目前 约 有 130 颗 CubeSat 3U (10x10x 
30cm) 卫星 *。 目 前 已 经 搭载 了 7 第 3 代 传 感 器 PSB. 
SD, 幅 宽 为 32.5x19.6 km, 能 采集 8 波段 的 多 光谱 影 
像 , 空 间 分 辩 率 为 3 m。 这 些 卫星 大 多 位 于 太阳 同 


1 研究 区 概况 图 
Fig. 1 Overview map of the study area 
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步 轨道 ,所 有 卫星 组 成 的 卫星 星座 PlanetScope 拥 有 
天 对 地 球 成 像 的 能 力 。 本 研究 使 用 的 Planet- 
Scope 波段 参数 见 表 1, 使 用 蓝 波段 (B2) 、 绿 波段 工 
(B4) 、 红 波段 (B86) 、 红 边 波段 (B7) 和 近 红 外 波段 
(B8) 进 行 了 后 续 人 研究 。 本 文 使 用 的 是 PlanetScope 
Level-3B 数据 ,过 境 时 间 为 2021 年 7 月 24 日 ,Level- 
3B 产 品 使 用 6S 辐射 传输 模型 和 MODIS 数据 作为 畏 
助 对 PlatnetScope 获 取 的 遥感 影像 进行 大 气 校正 ,并 
经 过 了 传感器 校正 和 正 射 校正 。 


表 1 PlanetScope 影 像 光谱 波段 参数 


Tab.1 Spectral band parameters of PlanetScope images 


波段 波段 名 波长 /nm 
蓝 波 段 B2 465-515 
绿 波段 B3 513~549 
绿 波段 工 B4 547~583 
黄 波 段 B5 600~620 
红 波段 B6 650~680 
红 边 波段 B7 697-713 
近 红 外 波段 B8 845-885 


将 光谱 指数 作为 预测 土壤 特性 模型 的 输入 变 
量 时 ,能 有 效 降 低 光 谱 反 射 率 误 差 ”。 通 过 波段 运 


算得 到 的 盐 度 指数 与 植被 指数 被 证 明 能 准确 反 演 
土壤 盐分 "“。 本 研究 使 用 PlanetScope 原始 波段 数 
据 计算 光谱 指数 ,计算 公式 如 表 2 所 示 。 
14 mRMR 特征 变量 选择 算法 

mRMR 是 一 种 过 滤 型 特征 选择 算法 ,其 算法 原 
理 简单 且 对 数据 的 要 求 低 ,在 土壤 盐 度 监测 领域 还 
没有 相关 的 研究 探讨 该 算法 的 应 用 。 因 此 ,本 研究 
结合 mRMR 方法 构建 了 土壤 盐 度 预测 模型 ,为 土壤 
盐 度 监测 领域 提供 了 参考 。mRMR 考虑 了 特征 之 
间 相 似 性 造成 的 数据 宛 余 问 题 ,通过 计算 特征 与 离 
散 化 电导 率 数据 类 别 之 间 的 互信 息 ,特征 与 特征 之 
间 的 互信 息 ,来 度量 它们 之 间 的 相关 性 。 然 后 , 保 
留 与 离散 化 电导 率 数 据 类 别 之 间 互 信息 值 较 高 的 
特征 ,同时 去 除 特征 之 间 互 信息 值 较 高 导致 数据 宛 
余 的 特征 六 。 任 意 2 个 特征 x 和 y 的 互信 息 计算 公 
式 为 : 


I(x,y)- | P, )iog dy (1) 


式 中 :7 为 特征 x 和 7 的 互信 息 ; P(r 为 特征 x 和 
y 的 联合 概率 密度 函数 ;P() 为 特征 x 的 边缘 概率 密 
度 函 数 ;P(y) 为 特征 y 的 边缘 概率 密度 子 数 。 特 征 与 
离散 化 电导 率 数据 类 别 之 间 的 互信 息 计 算 公 式 为 : 


表 2 基于 PlanetScope 影 像 光 谱 指数 计算 


Tab. 2 Spectral indices calculation based on PlanetScope images 


光谱 指数 名 称 公式 参考 
土壤 调节 植被 指数 SAVI SAVI = 1.5 x(NIR — Red)/(NIR + Red + 0.5) 37] 
比值 植被 指数 RVI RVI - NIR/ Red 38] 
增强 植被 指数 EVI EVI=2.5x(NIR— Red)/(NIR + 6 x Red - 7.5 X Blue + 1) 39] 

归 一 化 差 值 植被 指数 NDVI NDVI=(NIR - Red)/(NIR + Red) 40] 
归 一 化 差异 盐 度 指数 NDSI NDSI «(Red - NIR)/(Red + NIR) 41] 
亮度 指数 BI BI 2 JRed + NIR? 41] 

盐 度 指数 SI SI= Blue x Red 41] 

盐 度 指数 SH SII = Green x Red 42] 
盐 度 指数 SI2 SI2= JGreen’ + Red’ + NIR? 42] 
盐 度 指数 SI3 SI3 = JGreem + Red? 42] 

盐 度 指 数 S1 S1= Blue/Red 43] 
盐 度 指数 S2 S2 =(Blue-Red/Blue+ Red) 43] 

盐 度 指数 S3 S3 =(Greenx Red)/Blue 43] 
盐 度 指数 S4 S4 = JBluexRed 43] 
盐 度 指数 S5 S5 =(Bluex Red)/Green 43] 
盐 度 指数 S6 S6 =(Red x NIR)/Green 43] 


注 :Red 为 红 波段 ;Blue 为 蓝 波段 ;Green 为 绿 波段 ;NIR 为 近 红 外 波段 。 
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IXS,p)- is Y (xp) 究 。Bagging 在 训练 集中 进行 Bootstrap HFE , BI m 
MZ (2) — 个 样本 中 进行 随机 的 \ 有 放 回 的 采样 ,得 到 n 个 与 原 
max D(S.p) 


式 中 :DG 四 为 特征 与 类 别 之 间 的 平均 互信 息 ;3 为 
变量 集合 ;ye 起 为 特征 :与 类 别 p 之 间 的 互信 息 ;w 
为 特征 i;p 为 离散 化 电导 率 数据 类 别 。 特 征 与 特征 
之 间 的 互信 息 计算 公式 为 : 

» - s? Y Tr) 

min A(S) 
式 中 :R(S) 为 特征 与 特征 间 的 平均 互信 息 ;1(xsw) 为 特 
征 i 与 特征 j 之 间 的 互信 息 ;% 为 特征 i;% 为 特征 j。 
结合 公式 (2) 与 公式 (3) 得 到 最 大 相关 最 小 元 余 的 
判断 条 件 商 标准 : 
Visa - DIR 
max ¢(D,R) 

RP so (DR Jt UB nlii D RE SRZ 
的 互信 息 ;R 为 特征 与 特征 间 的 互信 息 。mRMR 算 
法 通过 Matlab 2019b 运行 fcmrmr 函数 实现 ,fscmrmr 
函数 自动 将 电导 率 数据 离散 化 处 理 。 
1.5 建 模 策 略 及 评价 
1.5.1 Bagging 回归 算法 装 袋 算法 (Bootstrap aggre- 
gating, Bagging) 是 Breiman“ 提 出 的 一 种 并 行 集成 
学 习 方 法 ,将 多 个 单一 的 预测 需 聚 合 在 一 起 ,对 结 
果 进 行 平均 处 理 , 进行 分 类 或 者 回归 预测 相关 的 研 


(3) 


(4) 


MLFIAK — 基 模 型 1 预测 1 


台 训 练 集 样本 个 数 相 等 的 采样 集 。 单 一 预测 需 基 
于 这 些 采 样 集 进行 独立 的 训练 ,在 进行 回归 分 析 
时 ,将 n 个 预测 器 的 均值 作为 综合 预测 结果 ,过 程 如 
图 2 所 示 。Bagging 回归 算法 有 效 降 低 了 训练 集 的 
方差 ,避免 预测 结果 出 现 过 拟 合 。 通 过 Python 3.9 
导入 机 需 学 习 库 sklearm 构建 Bagging 回 归 模 型 。 
1.5.2 模型 构建 及 评价 并 非 所 有 特征 都 能 提升 
Bagging 回归 算法 预测 性 能 ,不 相关 特征 和 宛 余 特 征 
往往 会 降低 模型 预测 的 精度 和 效率 。 为 了 降低 负 
面 影响 ,探索 PlanetScope 多 光谱 数据 预测 土壤 盐分 
的 潜力 ,所 有 特征 构成 集合 I ,mRMR 算法 选择 前 
30% 的 特征 构成 集合 I ,基于 集合 I 和 集合 卫 ,开发 
了 两 种 模型 Model-1 和 Model-T ( 表 3)。 将 采集 的 
84 个 土壤 样品 进行 随机 划分 (在 Python 3.9 中 
test_size=0.3 用 于 划分 数据 集 )。 所 有 样品 的 70%， 
约 58 个 土 样 被 划分 为 训练 集 , 用 于 训练 模型 ;所 有 
样品 的 30%, 约 26 个 土 样 被 划分 为 验证 集 , 用 于 模 
型 的 验证 。 为 了 确保 模型 的 稳健 性 ,提高 模型 的 适 
用 性 ,降低 随机 划分 样本 对 模型 预测 性 能 的 影响 ， 
本 研究 随机 划分 土壤 样本 3 次 ,获得 3 组 训练 集 和 
验证 集 (A.B 和 C)。 为 了 确保 模型 的 可 比 性 ,每 次 
划分 训练 集 和 验证 集 时 ,对 于 土壤 样本 与 集合 1 和 


| 子 训练 集 2 |—9| ” 基 模 型 2 预测 2 


训练 集 


子 训练 集 n 


子 抽样 3  1—9| 子 训练 集 3 — 基 模 型 3 预测 3 


-一 一 一 一 一 A 
s a Gael 


验证 集 


| 


图 2 Bagging 


回归 算法 预测 过 程 


Fig.2 Prediction process of Bagging regression algorithm 


表 3 PlanetScope 数 据 建 模 策略 细节 
Tab. 3 Details of the PlanetScope modeling strategies 


特征 变量 类 型 


模型 模型 方法 ”特征 选择 
Model-I Bagging d 
81,82,83,84,85,86) 
Model-ll Bagging mRMR 


波段 反射 率 数 据 (B2、B4、B6、B7、B8)+ 光 谱 指 数 (SAVI、RVI、EVI、NDVI、NDSI、BI.、SI、SIl ,SI2,SI3, 


波段 反射 率 数据 (B4)+ 光 谱 指 数 (EVI.S6.、S1 .SI2、NDSI.BI) 
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集合 了 构成 的 数据 集 , 划 分 模式 保持 一 致 。 

为 了 评估 Bagging 回归 模型 预测 EC 的 能 力 ,使 
用 决定 系数 ( 尼 ) . 均 方 根 误差 (RMSE) 和 四 分 位 数 的 
相对 预测 误差 (RPIQ)3 个 指标 进行 衡量 。 模 型 预测 
的 准确 性 与 记 .RPIQ 成 正比 ,与 RMSE 成 反比 , 即 局 
与 RPIQ 越 大 ,RMSE 越 小 模型 预测 性 能 越 好 ,计算 
公式 如 下 : 


(5) 


(6) 


(7) 


IQR=03-01 (8) 
式 中 ;$5, 为 预测 值 ;p; 为 实际 测量 值 ; 5 为 测量 值 的 
平均 值 ;n 为 采样 点 数量 ;IQR 为 验证 集中 测量 值 的 
四 分 位 间距 ;RMSE; 为 验证 集 的 均 方 根 误差 ;03 为 
上 四 分 位 值 ;01 为 下 四 分 位 值 。 


2 结果 与 分 析 


24 土壤 电导 率 描述 性 统计 
EC 作为 土壤 盐分 的 参考 ,EC 的 描述 性 统计 结 
果 如 图 3 和 表 4 所 示 。 整 个 EC 的 范围 为 0.35~ 


国 加 训练 集 A 加 训练 集 B 四 训练 集 C 国正 全 集 
OMERA [验证 集 B 验证 集 C 
140r — pš o- 均值 。 样本 


y y D Q 
数据 集 


注 :EC 为 电导 率 。 下 同 。 
图 3 电导 率 统计 图 


Fig.3 Statistical diagram of electrical conductivity 


表 4 GEEA BICHER EC 描述 性 统计 
Tab. 4 Descriptive statistics of datasets A, B, and C 


compared to the entire dataset EC 


数据 集 平均 值 /4$.m” 标准 差 /d$.m” 变异 系数 

训练 集 A 39.47 33.55 0.85 
验证 集 A 34.86 32.60 0.93 
训练 集 B 39.19 34.45 0.88 
验证 集 B 35.51 30.48 0.86 
训练 集 C 37.88 33.82 0.89 
验证 集 C 38.42 32.21 0.84 

全 集 38.04 33.13 0.87 


124.90 dS .m ,平均 值 为 38.04 dS- m~’, PREE N 
33.13 dS.m', 变异 系数 为 0.87, 为 中 等 变异 。EC 描 
述 性 统计 结果 显示 ,土壤 样本 随机 划分 3 次 得 到 和 A、 
B 和 C 数 据 集 ,它们 的 训练 集 和 验证 集 分 布 与 整个 
数据 集 的 分 布 相 似 。 因 此 ,训练 集 和 验证 集 的 样本 
可 以 代替 整个 数据 集 的 样本 进行 建 模 和 验证 。 
2.2 特征 变量 优选 

集合 I 由 5 个 波段 反射 率 数 据 B2、.B4 B6, B74 
B8,4 个 植被 光谱 指数 RVI、EVI、SAVI 和 NDVI,12 
个 土壤 盐 度 光谱 指数 S1、S2、S3、S4、S5、S6.、SI、SIl、 
SI2 SI3 .BI 和 NDSI 组 成 。 从 集合 [中 选择 最 能 反映 
土壤 盐 涡 化 发 生发 展 状 况 的 特征 变量 构建 土壤 盐 
渍 化 预测 模型 ,能 有 效 提升 模型 效率 和 预测 精度 。 
目前 mRMR 特征 选择 方法 还 未 广泛 用 于 土壤 盐 污 
化 监测 领域 ,本 研究 使 用 Matlab 2019b 统计 和 机 器 
学 习 工 具 箱 执行 人 cmrmr 函数 ,通过 mRMR 算法 计 
算 互信 息 , 并 依据 互信 息 商 标准 保留 最 相关 特征 ， 
过 滤 宛 余 特 征 。 经 过 互信 息 商标 准 选择 前 30% 特 
征 作为 集合 了 ,包括 1 个 波段 反射 率 数据 B4,1 个 植 
被 光谱 指数 EVI,4 个 土壤 盐 度 光谱 指数 S6、S1、SI2 
和 NDSI, 图 4 为 波段 反射 率 与 光谱 指数 特征 重 
要 性 。 
2.3 预测 模型 和 性 能 评估 

为 了 提高 模型 的 预测 精度 和 泛 化 程度 ,使 用 基 
于 步 长 为 1 的 网 格 搜索 ,对 模型 影响 较 大 的 基础 学 
习 髓 数量 进行 了 1000 次 连续 和 欠 代 调整 ,根据 模型 尼 
值 确定 模型 参数 最 优 值 。 在 进行 参数 优化 时 ,其余 
参数 为 默认 值 。 图 5 显示 了 模型 尼 与 模型 参数 在 连 
续 迭 代 过 程 中 的 变化 关系 , 当 屁 最 大 时 确定 基础 学 
习 器 数量 数值 ,重复 此 过 程 直到 确定 参数 最 优 值 ， 
Model-I 模型 基础 学 习 需 数量 为 4,Model-I 模型 基 
础 学 习 器 数量 为 9。 
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注 :SAVI 为 土壤 调节 植被 指数 ; RVI 为 比值 植被 指数 ;EVI 为 增强 
植被 指数 ;NDVI 为 归 一 化 差 值 植被 指数 ;NDSI 为 归 一 化 差异 盐 度 
指数 ;BI 为 亮度 指数 ;SI SI1 ,SI2 .SI3 81,82 ,83,84,85 
和 S6 为 盐 度 指数 。 
图 4 PlanetScope 波 段 反 射 率 和 光谱 指数 的 特征 重要 性 


Fig.4 Importance of PlanetScope band reflectance 


and spectral indices 
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基于 优化 后 的 参数 ,将 集合 I 和 和 集合 卫 随 机 划 
分 的 训练 集 作为 自 变 量 输入 Bagging 模 型 ,构建 EC 
预测 模型 Model-I 和 Model-I ,利用 随机 划分 的 验证 
数据 集 对 模型 的 预测 性 能 进行 验证 。 模 型 Model-I 
的 尼 与 RPIQ 均 高 于 Model-I ,RMSE 均 低 于 Model-I , 
这 说 明 mRMR 算法 可 以 选择 更 重要 的 特征 代替 全 
部 特征 ,降低 数据 的 复杂 性 和 宛 余 性 ,建立 更 准确 
和 更 高 效 的 模型 ( 表 5)。 在 随机 划分 的 三 组 数据 集 
ABI C , Model-II 5j Model-I 4H EU , ?2E35 3 
2796 , RPIQ 平均 提高 了 20% , RMSE 平均 降低 了 
16%。mRMR 在 选择 特征 数量 较 少 的 条 件 下 ,能 
留 最 相关 最 小 宛 余 特征 ,能 最 大 化 特征 包含 的 有 用 
信息 。 在 集合 I 中 光谱 指数 EVI 和 S6 重 要 性 最 高 ， 
表明 EVI 对 土壤 盐 渍 化 引起 的 植被 状况 与 光谱 反射 
率 变化 敏感 ,S6 对 盐 渍 土 表 面 状况 与 光谱 反射 率 变 
化 敏感 。 光 谱 指 数 EVI 与 S6 在 土壤 盐 渍 化 监测 中 
具有 巨大 潜力 。 


ü 


0.80 r (a) Model-I 的 优化 结果 0.80- (b) Model-I 的 优化 结果 
0.75 0.75 - 
0.70 0.70 +- 
0.65 0.65 - 
R 0.60 & 0.60 - 
0.55 0.55 L 
0.50 0.50 - 
0.45 dust. 
940p LU o 040 
0 200 400 600 800 1000 1200 0 200 400  Á 600 $800 1000 1200 
迭代 次 数 迭代 次 数 
注 : 尼 为 决定 系数 。 下 同 。 
图 5 模型 参数 调整 
Fig. 5 Model parameters adjustment 
ms 模型 预测 性 能 比较 
Tab. 5 Comparison of model prediction performance 
数据 集 模型 uii. i RPIQ 
R RMSE/dS -+ m” R RMSE/dS- m” 
A Model-I 0.79 15.26 0.67 18.41 3.34 
Model-II 0.90 10.27 0.76 15.78 3.93 
B Model-[ 0.79 15.82 0.41 22.95 227 
Model-II 0.87 12.48 0.66 17.45 3.00 
C Model-[ 0.76 16.43 0.49 22.58 2.42 
Model-II 0.87 11.93 0.57 20.75 2.71 
平均 Model-[ 0.78 15.84 0.52 21.32 2.68 
Model-II 0.88 11.56 0.66 18.00 3.21 


主 : 尺 为 决定 系数 ;RMSE 为 均 方 根 误差 ;RPIQ 为 四 分 位 数 的 相对 预测 误差 。 
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2.4 土壤 盐 度 空间 分 布 

基于 模型 Model-[ 和 Model-I ,使 用 随机 划分 的 
数据 集 A、B 和 C 绘 制 了 3 张 高 分 辨 率 的 土壤 EC 分 
布 图 ,对 模型 Model-[ 和 Model-I 的 3 张 预测 图 进行 
平均 ,得 到 网 6a 和 图 6b。 图 6a 和 图 6b 产 生 了 相似 
的 土壤 EC 空间 分 布 ,Model-T 模型 产生 了 最 理想 的 
结果 ,预测 的 土壤 EC 范围 (0.70~100.83 dS.m ') 更 符合 
实际 测量 的 土壤 EC 值 。 图 6b 相 较 图 6a, 更 准确 的 
显示 了 农田 区 域 与 荒漠 区 域 土壤 EC 分 布 的 细微 变 
化 ,更 符合 该 区 域 土 壤 盐 度 分 布 。 受 地 形 影 响 , 土 
坏 盐 度 高 值 主 要 分 布 在 荒漠 区 域 ; 受 地 形 和 灌溉 的 
影响 ,土壤 盐 度 低 值 主 要 分 布 在 农田 区 域 。 


3 讨论 


PlanetScope 数据 结合 mRMR 特征 选择 方法 首 
次 用 于 绘制 高 空间 分 辩 率 的 土壤 EC 图 , ES i 
化 监测 领域 缺乏 类 似 的 研究 作为 参考 ,因此 借鉴 了 
马 慧 右 等 ”和 Jing 等 “的 研究 ,通过 比较 模型 
Model-I fill Model-II PHW RE CK 5) ,得 出 与 其 
相似 的 研究 结果 ,mRMR 方法 能 优化 建 模 特征 并 改 
善 模型 预测 性 能 。 基 于 PlanetScope 数据 建立 的 模 
型 Model-I 和 Model- 荆 获取 的 预测 值 如 图 7 所 示 ， 
Model-T 获 取 的 预测 值 与 实测 值 更 接近 1:1 拟 合 线 ， 
表明 Model-T 能 获取 更 准确 的 预测 结果 。 基 于 Plan- 
etScope 影 像 低 于 $m 的 空间 分 辩 率 ,能 避免 如 Scud- 
iero 等 “使 用 Landsat7 多 光谱 影像 进行 土壤 盐 度 评 
估 时 ,结果 难以 表征 土壤 盐 度 的 空间 变异 性 等 情 
况 。 也 能 避免 使 用 MODIS 系列 Landsat 系列 或 Sen- 
tinel 系列 数据 进行 土壤 盐 度 监测 时 ,预测 结果 难以 


(a) 基于 Model-I 的 EC 空间 分 布 


ES 


人 


zm 


(b) 基于 Model-I 的 EC 空间 分 布 


具备 指导 意义 等 情况 ”。 基 于 此 ,使 用 PlanetScope 
数据 结合 mRMR 特征 选择 方法 建立 的 模型 Model- M 
对 研究 区 土壤 EC 进行 预测 并 绘制 成 图 ,避免 了 由 
于 空间 分 辩 率 不 足 对 土壤 盐 度 监测 和 绘制 带 来 的 
限制 。 

本 研究 将 光谱 指数 作为 自 变 量 输入 Bagging 模 
型 ,改善 模型 的 预测 性 能 ,提高 模型 定量 描述 土壤 
盐 度 的 能 力 , 丁 建 丽 等 “与 曹 雷 等 ”的 研究 结果 证 
实 这 是 可 行 的 。 并 非 所 有 光谱 指数 都 能 用 于 优化 
模型 ,其 中 不 相关 和 和 宛 余 的 特征 变量 反而 会 起 反 作 
用 。 本 文通 过 mRMR 方法 对 所 有 特征 变量 进行 了 
重要 性 排序 ,将 排序 后 的 特征 变量 作为 自 变量 输入 
模型 Model- 中 。 当 选择 的 特征 变量 数 小 于 309% 
时 ,模型 的 预测 精度 较 低 ; 当选 择 的 特征 变量 数 超 
过 30% 时 ,模型 的 预测 性 能 略 低 基本 保持 稳定 。 

重要 性 排序 中 FEVI 重 要 性 最 高 , 远 远 高 于 其 他 
波段 与 光谱 指数 (图 4) ,对 土壤 盐 涡 化 具备 更 多 的 
解释 能 力 。 这 可 能 是 因为 植物 叶 面 积 和 冠 层 结构 
的 发 育 受 到 土壤 盐 汗 化 抑制 更 明显 ,而 EVI 对 于 植 
物 冠 层 结构 变化 更 敏感 , 故 EVI 在 建 模 时 更 重 
Zi. Rao 等 ”的 研究 表明 , 盐 渍 十 具备 更 高 的 光谱 
反射 率 , 且 在 可 见 光 和 近 红 外 波长 范围 , 盐 渍 化 程 
度 与 光谱 反射 率 成 正比 ,在 盐 渍 化 程度 较 高 的 农田 
苇 漠 接壤 处 与 荒漠 区 域 可 以 直接 利用 PlanetScope 
原始 波段 反射 率 对 土壤 盐 度 进行 建 模 分 析 。mRMR 
选择 的 特征 集合 正中 仅 包 含 原 始 波段 B4, 故 B4 相 较 
于 其 他 原始 波段 ,对 于 土壤 盐 渍 化 更 敏感 ,更 适合 
用 于 描述 土壤 盐 度 信 息 。 在 本 研究 中 除 植被 覆盖 
度 较 高 的 农田 以 外 , 均 能 使 用 对 土壤 盐 度 敏 感 的 盐 
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图 6 土壤 电导 率 空间 分 布 


Fig.6 Spatial distributions of soil electrical conductivity 
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(a) 基于 数据 集 A 的 Model-I 预 测 (b) 基于 数据 集 B 的 Model-I 预 测 (c) 基于 数据 集 C 的 Model-I 预 测 
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iE: RMSE 为 均 方 根 误差 ;RPIQ 为 四 分 位 数 的 相对 预测 误差 。 
图 7 模型 ModelLI 和 Model-I 预 测 性 能 比较 


Fig.7 Comparison of prediction performance of model Model-[ and Model-II 


度 指 数 定量 描述 土壤 盐 涡 化 情况 ,mRMR 选择 的 前 
30% 特 征 变量 中 ,S6、B4、S1 和 NDSI 均 为 盐 度 指数 ， 
故 盐 度 指数 对 于 建立 土壤 盐 度 预测 模型 必 不 可 
少 。 汉 娟 等 5 Mettemicht 等 中 和 丁 建 丽 等 3 的 研 
究 结 果 证 明了 这 些 观点 是 合理 的 。 

土壤 样本 的 随机 划分 和 模型 参数 的 设 定 使 模 
型 的 误差 出 现 累积 。 随 机 划分 的 训练 集 和 验证 集 
与 整个 数据 集 的 分 布 相似 , 可 以 代替 整个 数据 集 的 
样本 进行 建 模 和 验证 。 为 避免 样本 随机 划分 为 模 
型 带 来 的 不 确定 性 ,将 随机 划分 3 次 得 到 的 数据 集 
A.B 和 C 用 于 训练 模型 和 验证 模型 ,最 终 对 评价 指 
标 取 均值 来 评估 模型 ,增加 了 模型 预测 的 稳健 性 ， 
Jing 等 外 的 研究 结果 证 明 这 是 有 效 的 。 虽 然 多 次 随 
机 划分 样本 能 有 效 降 低 模 型 预测 的 不 确定 性 ,但 是 
本 研究 并 没有 确定 最 有 效 的 随机 划分 次 数 。 在 设 
定 模 型 参数 时 ,为 了 降低 模型 参数 对 预测 性 能 的 影 
响 ,基于 步 长 的 网 格 搜索 ,对 模型 影响 较 大 的 基础 
学 习 带 数量 进行 了 连续 迭代 调整 ,结果 如 图 5 所 
示 。 基 于 步 长 的 网 格 搜索 是 否 是 参数 优化 的 最 佳 
方式 ,这 也 存在 不 确定 性 。 除 上 述 两 点 ,人 类 活动 
也 增加 了 模型 预测 的 不 确定 性 ,图 6a 和 图 6b 左 上 角 


区 域 显示 ,通过 Model-[ 和 Model-I 预测 的 农田 区 域 
土壤 盐 度 被 高 估 ,这 与 实际 情况 不 符 。 通 过 野外 调 
查验 证 ,该 区 域 为 了 达到 节 水 增产 的 目的 ,推广 了 
膜 下 滴灌 技术 。 履 膜 区 域 光 谱 反 射 率 不 同 于 裸 土 
区 域 , 整 体 偏 高 ,由 于 多 光谱 影像 的 光谱 分 辩 率 不 
AE. ,使 得 盐碱地 与 覆 膜 地 难以 进行 区 分 ,在 未 来 的 
人 研究 中 将 尝试 利用 高 光谱 技术 来 消除 农田 区 域 覆 
膜 给 土壤 盐 度 监测 带 来 的 影响 。 此 外 ,本 研究 仅 考 
IŠ T Bagging 回归 模型 , 马 国 林 等 “的 研究 结果 表 
明 , 仅 通 过 一 个 模型 难以 取得 最 佳 预测 效果 。 
此 ,在 未 来 的 研究 中 将 构建 多 种 模型 对 土壤 盐 度 进 
行 预测 ,并 对 比 模型 性 能 选择 最 适合 本 研究 数据 集 
的 模型 。 


4 结论 


本 人 研究 基于 PlanetScope 影 像 原始 波段 数据 、 植 
被 光谱 指数 和 土壤 盐 度 指数 共 21 个 特征 变量 ,使 用 
mRMR 方法 筛选 前 30% 特 征 变量 并 输入 Bagging 回 
归 模 型 ,建立 Model-I 和 Model-I 土壤 盐 度 预测 模 
型 ,成 功 绘制 了 研究 区 3 m 2: SEXE I] EC 空间 分 布 
图 ,为 准确 监测 随时 空 动态 快速 变化 的 土壤 盐 度 信 
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息 提供 参考 。 得 出 以 下 结论 : 

(1) 将 PlanetScope 影 像 原始 波段 数据 .植被 光 
谱 指 数 和 土壤 盐 度 指数 共 21 个 特征 变量 作为 自 变 
量 建立 的 Model-I 模型 ,基本 可 以 实现 对 EC 的 预测 
(验证 集 平均 尺 =0.52, 平 均 RMSE=21.32 dS.m', 平 
均 RPIQ=2.68)。 

(2) 将 mRMR 筛选 后 的 特征 变量 作为 自 变量 
立 的 Model-T 模 型 ,降低 了 预测 EC 的 特征 维度 , 提 
高 了 模型 的 预测 效率 和 准确 性 (验证 集 平均 R= 
0.66 ,平均 RMSE=18.00 dS.m', 平 均 RPIQ=3.21)。 
Model-I 与 Model-I 相 比 ,EVI、S6、B4、S1、SI2 和 
NDSI 对 农田 荒漠 交错 带 士 壤 盐 涡 化 监测 具有 重要 
作用 。 
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Digital mapping of soil salinization in a typical oasis based 
on PlanetScope images 


LI Ke", DING Jianli", HAN Lijing", GE Xiangyu", GU Yongsheng", 
ZHOU Qian", LYU Yangxia 
(1. Xinjiang Common University Key Lab of Smart City and Environment Stimulation, College of Geography and Remote Sensing 
Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 2. Key Laboratory of Oasis Ecology under Ministry 


of Education, Xinjiang University, Urumqi 830046, Xinjiang, China; 3. Groundwater and Salt Monitoring Station, 
Ogan River Basin Management Office, Aksu 842000, Xinjiang, China) 


Abstract: High-resolution soil salinity maps are urgently needed in arid and semi-arid regions to visualize the 
subtle spatial variations in salinity distribution. These maps are crucial for guiding the development of land re- 
source management policies and water resource management policies in salt-affected and potentially salt-affected 
areas, aiming to prevent further soil degradation and ensure sustainable agricultural economic development and 
food security. Based on PlanetScope imagery, vegetation spectral indices and soil salinity indices were extracted, 
resulting in a total of 21 variables. These variables were used as input for the Bagging algorithm to construct a 
soil salinity prediction model, referred to as Model-I . The max-relevance and min-redundancy (mRMR) method 
was employed to select relevant feature variables, which were then inputted into the Bagging algorithm to build a 
soil salinity prediction model, referred to as Model-II . Field sampling data were used to assist in model building 
and validation. Model-I and Model-Il were evaluated using model evaluation metrics. The results indicate that 
the prediction performance of Model-I is better than that of Moedl-I (mean &^—0.66, mean RMSE-18.00 dS* m`, 
mean PRIQ-3.21 for the validation set), and that mRMR effectively reduces the multidimensional feature redun- 
dancy. PlanetScope images combined with the mRMR method successfully mapped high-resolution soil salinity, 
which provided more detailed information on the spatial distribution of soil salinity, and the results of the study 
promoted the use of PlanetScope data to monitor soil salinity information. 


Key words: PlanetScope; mRMR; Bagging; soil salinization; digital mapping techniques 


